PUDLEY  KNOX  LIBRARY 
NAVAL  POSTGT  PE  SCHOOL 

£8-6008 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


A 

COMPARISON  OF  SATELLITE-DERIVED  OCEAN 

VELOCITIES 

WITH  OBSERVATIONS  IN  THE  CALIFORNIA 

COASTAL  REGION 

by 

John  F.  O'Hara 

March  1987 

Thesis 

Advisor         David  C.  Sm 

ith 

IV 

Approved  for  public  release;  distribution  is  unlimited. 


T 233627 


c    'm   C    PAGE 


REPORT  DOCUMENTATION  PAGE 


<s    REPORT   SECuRiTy  CLASSIFICATION 

UNCLASSIFIED 


'D     RESTRICTIVE    MARKINGS 


Ja    SECUR'TY  CLASSIFICATION  AUTHORITY 


.D    OEl^aSS  c  CAT  ON     DOWNGRADING   SCHEDULE 


3      DISTRIBUTION/ AVAILABILITY   Of    REPORT 

Approved  for  public  release; 
distribution  is  unlimited. 


PERFORM  \G   ORGANIZATION   REPORT   NuMBER(S) 


S    MONiTORiNG   ORGANIZATION   REPORT   NJV3ER(Sj 


ea     SAME   OF   PERFORMING  ORGANIZATION 


Naval  Postgraduate  School 


60    OFFiCE    SyMBOL 
(if  jppiiable) 

Code    63 


?a    NAME   OF   MONiTORiNG   ORGANIZATION 

Naval  Postgraduate  School 


be    ADDRESS    C<f>    State    *nd  /iP Code) 

Monterey,  California  93943-5000 


7b     ADDRESS  (Ory.   State    try)  ZIP  Code) 
Monterey,    California    93943-5000 


SAVE   0;    ^NDiNG- SPONSORING 
ORGAVZAT  ON 


8b    OFFiCE  SYMBOL 
(it  applicable) 


9    PROCUREMENT   INSTRUMENT   lDEN".f  XATlON    NUMBER 


8c    ADDRESS  iGry    State  and  HP  Code) 


10    SOURCE    OF   FUNDING   NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


rAS-C 

NO 


WORK    JNIT 
ACCESS;ON   NO 


'.£   tnc'uOe   Security  Classification) 
A   COMPARISON   OF    SATELLITE-DERIVED   OCEAN   VELOCITIES   WITH   OBSERVATIONS    IN   THE 

CALIFORNIA  COASTAL  REGION 


•  A.  A, 


:R«S: 


O'Hara,  John,  F. 


y  r.  tc  at  pOR" 

Master's   Thesis 


"3d    "ME  c  o  «  e  r  e  d 

tOQM  T0 


fl    DATE  OF   REPORT    (Year   Month   Day) 

1987   March 


'S    -AGE    CD. 
42 


n:-.iV:','A9'     SO'ATQN 


:OSA'    CODES 


SuBGROL'P 


'8    SuBjECT    TERMS  'Continue  on  reverse   ''  neceua^y   and  identify   by   bicck    number) 

Remote  Sensing,  Ocean  Currents,  California  Current 
System,  California  Coastal  Region 


"3    ~d'>rRAC~     Continue  on  reverie   it  necessary   tnd  identity  by   block   number) 

Satellite-derived  ocean  surface  velocity  vectors  in  the  California 
Current  System  (CCS)  are  compared  with  in  situ  hydrographic  and  Doppler 
data.  The  in  situ  data  were  acquired  during  the  April  1981  phase  of  the 
Coastal  Ocean  Dynamics  Experiment  (CODE)  experiment.  In  general,  the 
satellite-derived  velocities  agreed  with  the  in  situ  data.  Due  to  the 
baroclinic  nature  of  the  study  region,  the  satellite  vectors  were  found 
to  be  representative  of  the  subsurface  geostrophic  flow.  Although  the 
number  and  concentration  of  the  satellite  vectors  was  small,  these 
vectors  were  capable  of  resolving  the  mesoscale  features  located  in  the 
study  region.  Comparison  of  co-located  data  revealed  that  the  in  situ 
Doppler  velocity  measurements  were  approximately  1.  5  times  larger  than 
the  satellite-derived  velocity  vectors.  These  results  agree  with  a 
similar  study  by  Flament  et  al.  (1985)  which  was  conducted  in  an 
offshore  region. 


20     j  S'»'3^T  ON-  AVAILABILITY  OF   ABSTRACT 

S§  -NClASS'F'EDUNI'MITED       □  SAME  AS   RPT  O  DTiC   USERS 


21     ABSTRACT  SECURITY   CLASSIFICATION 

UNCLASSIFIED 


.a    SAME    OF   RESPONSIBLE     NDiViOUAL 

David   C.    Smith    IV 


220  TELEPHONE  (include  Are*  Code) 

(408)    646-3350 


22c    OFFICE    SYMBOL 

Code    68Si 


3D  FORM  1473.84MAR 


83  apr  edit'on  '▼■ay  be  used  untie«tiausted 
All  otner  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF   th>S  PAGE 

UNCLASSIFIED 


Approved  lor  public  release;  distribution  is  unlimited. 


A  Comparison  of  Satellite-Derived  Ocean  Velocities 
with  Observations  in  the  California  Coastal  Region 


by 


John  F.  O'Hara 

Lieutenant,  United  States  Navv 

B.S.,  United  States  Naval  Academy^  1980 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  METEOROLOGY  AND  OCEANOGRAPHY 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 

March  1987 


ABSTRACT 

Satellite-derived  ocean  surface  velocity  vectors  in  the  California  Current  System 
(CCS)  are  compared  with  in  situ  hydrographic  and  Doppler  data.  The  in  situ  data 
were  acquired  during  the  April  1981  phase  of  the  Coastal  Ocean  Dynamics  Experiment 
(CODE)  experiment.  In  general,  the  satellite-derived  velocities  agreed  with  the  in  situ 
data.  Due  to  the  baroclinic  nature  of  the  study  region,  the  satellite  vectors  were  found 
to  be  representative  of  the  subsurface  geostrophic  How.  Although  the  number  and 
concentration  of  the  satellite  vectors  was  small,  these  vectors  were  capable  of  resolving 
the  mesoscale  features  located  in  the  study  region.  Comparison  of  co-located  data 
revealed  that  the  in  situ  Doppler  velocity  measurements  were  approximately  1.5  times 
larger  than  the  satellite-derived  velocity  vectors.  These  results  agree  with  a  similar 
study  by  Flament  et  al.  (1985)  which  was  conducted  in  an  offshore  region. 
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I.  INTRODUCTION 

As  the  cost  of  oceanographic  hydrographic  surveying  (i.e.  ship  time,  equipment, 
personnel,  etc.)  continues  to  rise,  satellite  remote  sensing  of  the  ocean  environment  has 
become  increasingly  popular.  Although  remote  sensing  does  have  some  limitations 
(e.g.  cloud  coverage  in  Sea  Surface  Temperature  (SST)  measurements),  satellites  do 
allow  oceanographic  data  to  be  collected  over  very  large  areas  in  an  extremely  short 
amount  of  time.  This  synoptic  view  of  a  large  region  is  not  routinely  available  in 
oceanographic  research  other  than  by  satellite.  Satellite  sensing  also  allows  researchers 
and  ordinary  users  to  "resample"  their  particular  area  of  interest  sometimes  several 
times  a  day.  This  availability  of  data  has  provided,  and  will  continue  to  provide,  the 
impetus  for  research  into  determining  how  to  acquire  the  maximum  amount  oC 
information  from  the  image(s).  A  particular  technique  for  extracting  satellite 
information  is  addressed  in  this  thesis.  Of  importance  for  the  Navy  is  the  ability  to  use 
satellite-derived  data  in  ocean  model  initialization  for  prediction  purposes.  Significant 
advances  could  also  be  made  in  weather  forecasting  using  an  air-ocean  coupled  model 
initialized  with  such  satellite  data. 

This  research  centers  on  comparing  satellite  data  with  in  situ  hydrographic  and 
Doppler  data.  It  is  the  object  of  this  paper  to  verify  coastal  ocean  flow  vector  fields, 
that  have  been  generated  using  a  feature  tracking  technique  on  three  sequential 
infrared  (IR)  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  images,  in  the 
California  Current  System  (CCS).  How  well  these  vectors  represent  the  coastal  ocean 
flow  will  be  determined  by  comparing  them  with  in  situ  data  collected  during  the 
Coastal  Ocean  Dynamics  Experiment  (CODE)  experiment.  Figure  1.1  shows  the  study 
region. 

Chapter  II  presents  background  information  on  the  general  character  of  the  CCS 
and  previous  investigations  that  have  addressed  satellite-derived  flow  vectors.  The 
techniques  used  to  generate  the  velocity  vectors  and  the  associated  error  budget  are 
described  in  Chapter  III.  Chapter  IV  compares  satellite  results  obtained  in  the  study 
with  hydrographic  and  Doppler  data.  Conclusions  and  recommendations  are  presented 
in  Chapter  V. 
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Figure  1.1     Box  indicates  CODL  satellite  comparison  region. 
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II.  BACKGROUND 

A.       PHYSICAL  OCEANOGRAPHY  OF  THE  CALIFORNIA  CURRENT  SYSTEM 

Hickey  (1979)  describes  the  climatological  characteristics  of  the  currents  which 
comprise  the  CCS;  the  California  Current,  the  California  Undercurrent,  the  Davidson 
Current  and  the  Southern  California  Countercurrent.  The  California  Current  is  defined 
as  a  surface  equatorward  flow,  while  the  California  Undercurrent  is  defined  as  a 
subsurface  poleward  flow  over  the  slope.  The  Davidson  Current  is  defined  as  a 
poleward  flow  that  occurs  during  the  fall  and  winter  north  of  Point  Conception.  The 
Southern  California  Countercurrent  is  also  defined  as  a  poleward  flow.  However,  it 
occurs  south  of  Point  Conception  and  inshore  of  the  Channel  Islands  in  the  California 
Bight.  Since  the  research  in  this  paper  centers  on  the  central  California  coastal  region, 
the  current  of  particular  interest  is  the  California  Current. 
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Fieure  2.1     Seasonal  profiles  of  the  meridional  California  Current 
oil' northern  CaHtornia  (from  Ikeda  and  Emery,  1984). 

The  seasonal  variation  of  the  California  Current  is  shown  schematically  in  Figure 
2.1  (Ikeda  and  Emery.  19S4).  There  are  three  axes  of  strong  southward  flow.  The 
nearshore  flow  is  most  fully  developed  in  the  spring  and  early  summer  south  of  Cape 
Blanco.  The  offshore  flow  on  the  other  hand,  is  most  fully  developed  in  the  late 
summer  or  fall.  South  of  Cape  Mendocino,  a  second  offshore  region  of  strong 
southward  flow  is  present  throughout  the  year.  The  nearshore  flow  is  driven  by  the 
longshore  component  of  the  wind  stress,  while  the  two  offshore  flows  are  driven  by  the 
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curl  of  the  wind  stress  (Hickey,    1979).     Various  counterflows  occur  north  of  Point 

Conception  near  the  coast.    During  the  fall  and  winter  these  counterflows  displace  or 

in  some  cases  even  replace  the  nearshore  southward  current.    While  the  large  scale 

climatology  of  the  California  Current  can  be  depicted  quite  easily,  the  daily  variations 

in  its  mesoscale  conditions  are  much  more  complex. 

To  investigate  these  mesoscale  features,  a  field  program  (CODE)  was  designed. 

The   basic  design   of  the   survey  grid  (Huyer  and   Kosro.    1987)  was  influenced   by 

previous  studies  of  other  upwelling  regions  and  as  a  result,  the  alongshore  spacing 

between  sections  was  much  larger  (20-30  km)  than  the  onshore-offshore  spacing  (5-10 

km).    It  should  be  noted  that,  because  of  this  spacing,  dynamic  topography  could  not 

resolve  eddies  or  meanders  with  alongshore  scales  less  than  40-60  km.    Each  section 

extended  about  45  km  from  the  coast  (spanning  the  continental  shelf  and  upper  slope) 

and  each  of  the  sections  were  chosen  to  be  locally  perpendicular  to  the  shelf  isobaths. 

Conventional  CTD  casts  to  a  maximum  depth  of  1000  m  and  continuous  operation  of 

a  Doppler  acoustic  log  (DAL)  to  a  depth  of  150  m  (to  obtain  velocity  profiles)  were 

included  in  each  section.    The  geopotential  anomaly  (AD)  of  the  sea  surface  relative  to 

100  dbar  was  computed  for  each  survey.    For  the  few  stations  that  were  shallower  than 

100m.   ADn  was  calculated  directlv  and  AD  ,nn  was  calculated  bv  linear 

0  pmax  *  pmax  100 

extrapolation  from  the  next  offshore  pair  of  stations  to  obtain  AD0/100. 

Huyer  and  Kosro  (1987)  found  that  the  resulting  geopotential  anomaly  maps 
were  very  similar  to  maps  of  near-surface  salinity.  Lowest  values  of  geopotential 
anomaly  coincided  with  the  high-salinity  upwelled  water  near  shore,  and  the  highest 
values  were  associated  with  the  low-salinity  water  offshore.  The  geostrophic  flow 
appeared  to  be  along  isohalines,  with  the  low-salinity  water  on  the  right  and  generally 
equatorward  over  the  shelf  and  upper  slope.  During  the  April  1981  survey,  a  baroclinic 
jet  was  detected,  but  its  position  varied  widely  over  the  study  region.  The  strongest 
portion  of  this  jet  was  located  just  offshore  near  Point  Arena. 

Because  of  the  continuous  operation  of  the  DAL,  the  current  vectors  generated 
from  the  Doppler  acoustic  log  show  substantially  increased  mesoscale  variations  in  the 
flow  pattern  as  compared  to  the  geopotential  patterns.  Huyer  and  Kosro  (1987)  state 
that  maps  of  the  current  vectors  can  resolve  eddies  or  meanders  with  alongshore  scales 
of  at  least  20-30  km  (the  distance  between  sections).  Although  the  dynamic 
topography  of  the  April  1981  survey  indicates  a  weak  and  meandering  baroclinic  jet, 
the  associated  DAL  maps  do  not  show  clear  evidence  of  a  continuous  equatorward 
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current.  Rather,  during  the  April  19S1  survey  smaller-scale  eddy-like  features,  which 
penetrated  down  to  at  least  150m.  were  observed  over  the  slope  and  shelf.  A  more 
complete  discussion  of  the  April  19S1  survey  results  and  a  comparison  with  satellite 
data  are  presented  in  Chapter  IV. 

Modeling  studies  (Ikeda  and  Emery,  1984)  of  similar  mesoscale  circulation 
features  off  the  coast  of  Vancouver  Island  have  shown  that  the  growth  of  these  feature 
is  due  to  the  baroclinic  instability  associated  with  the  vertical  shear  between  the  surface 
current  and  the  undercurrent.  Ikeda  and  Emery  propose  that  the  mesoscale  features 
seen  in  satellite  imagery  off  the  California  coast  are  due  to  the  same  mechanism 
because  of  the  great  similarity  in  the  mean  current  structure  found  off  Vancouver 
Island  and  California.  After  modifying  the  nonlinear  numerical  model  for  the 
California  coast.  Ikeda  and  Emery  were  able  to  successfully  simulate  the  initial 
excitation,  growth,  and  cascade  of  energy  to  longer  wavelengths  of  these  features. 

Flament  et  al.  (1985)  studied  the  evolution  of  an  upwelling  filament  offshore 
along  the  coast  of  California  during  June  1981.  Velocity  vectors  were  generated  by 
tracking  the  cores  of  small  eddies  which  Flament  et  al.  assumed  to  be  advected  by  the 
large-scale  flow.  A  comparison  of  the  satellite-derived  velocity  vectors  against  in  situ 
DAL  current  vectors  revealed  that  the  Doppler  log  speeds  were  larger  by  a  factor  of 
approximately  1.5.  Flament  et  al.  proposed  that  either  the  speed  of  tracked  jet 
increased  during  the  hydrographic  survey,  or  there  was  a  systematic  bias  of  the 
satellite-derived  velocities. 

B.       REMOTELY  SENSED  OCEAN  DATA 

In  the  past,  AVHRR  data  has  been  used  to  intercompare  with  in  situ  sea  surface 
temperature  measurements  (Bernstein,  1982),  in  situ  dynamic  height  measurements 
(Rienecker  et  al.,  1985),  and  most  recently,  to  track  ice  movements  (Xinnis  et  al., 
19S6).  Feature  tracking  research  can  be  broken  into  several  different  approaches  based 
on  determining;  1)  how  to  track  features  that  appear  in  AVHRR  images;  and  2)  which 
features  to  track. 

Saunders  (1973)  used  an  airborne  IR  radiometer  to  estimate  the  SST  of  an  area 
90  km  X  90  km  off  the  New  England  coast.  This  area  was  remapped  over  a  two-week 
period  which  provided  enough  data  to  allow  both  time  derivatives  and  horizontal 
gradients  of  SST  to  be  computed.  Using  this  information,  Saunders  developed  a 
technique  to  compute  surface  flow.    However,  Saunders  had  to  make  the  following 
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assumptions:  that  the  mixed  layer  temperature  was  a  quasi-conservative  quantity,  and 
that  horizontal  diffusion  of  temperature  could  be  neglected  for  scales  of  5  km  and  time 
intervals  of  a  few  days.  The  technique  predicted  that  the  component  of  How  parallel  to 
the  mean  horizontal  temperature  gradient  was  solely  due  to  the  rate  of  change  of  the 
SST.  Unfortunately,  the  technique  could  not  predict  the  component  of  flow  normal  to 
the  temperature  gradient  unless  the  horizontal  flow  in  the  upper  ocean  was  assumed  to 
be  horizontally  non-divergent.  The  normal  velocity  component  was  then  computed  by 
determining  the  normal  displacement  of  an  isotherm  between  a  pair  of  flights.  Due  to 
a  lack  of  in  situ  surface  current  measurements,  Saunders  was  unable  to  determine 
exactly  how  accurately  this  method  represented  true  conditions.  Vet,  this  technique 
represents  a  start  in  trying  to  determine  whether  surface  flow  information  could  be 
obtained  from  a  remote  sensing  platform. 

Another  method  was  used  by  Kelly  (1983)  to  compare  satellite-derived  SST  maps 
with  in  situ  current  velocity  measurements  in  the  CCS.  Kelly's  method  required 
"decimation"  of  the  images  by  averaging  over  5x5  boxes  of  data.  The  "decimation"  of 
the  data  consisted  of  computing  EOF's  (Empirical  Orthogonal  Functions)  using  a 
computationally  faster  SVD  (Singular  Value  Decomposition)  method.  However,  the 
SVD  method  had  a  weakness  in  its  inability  to  handle  missing  data  and  required  that 
any  missing  data  (due  to  cloud  coverage)  be  filled  in.  These  data  gaps  were  filled  in 
using  a  linear  combination  of  data  in  space  and  time,  which  represented  the  best 
estimate  in  a  least-squares  sense.  Only  "small"  gaps  were  filled  by  spatial  interpolation 
due  to  Kelly's  requirement  that  a  missing  datum  could  only  be  replaced  if  the  mean 
square  error  of  the  estimate  was  less  than  30%  of  the  variance  at  that  location. 
Composite  images  were  created  using  partial  images  separated  by  12  or  24  hours  to  fill 
in  large  holes.  The  SST  EOF's  were  then  computed  using  these  interpolated  and 
composite  images.  Kelly  showed  that  the  EOF's  of  the  SST  field  revealed  that  the 
shapes  of  the  SST  mean  and  anomalies  were  dominated  by  local  effects  (coastal 
topography,  orographic  effects,  and  spatial  variations  in  the  wind).  The  EOF's  also 
revealed  that  the  intensity  of  the  anomalies  varied  with  the  strength  of  the  local  wind. 
Then  Kelly  showed  that  the  SST  patterns  were  the  same  as  those  of  the  near-surface 
velocities,  which  suggested  that  the  near-surface  current  field  may  have  also  been 
forced  by  the  local  variations  in  the  wind. 

Vastano  and  Reid  (1985)  used  an  interactive  algorithm  to  estimate  the  motion  of 
sub-mesoscale  features.    This  interactive  algorithm  required  the  operator  to  subjectively 
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generate  the  displacement  vectors  of  individual  features  in  successive  pairs  of  AVHRR 
images  at  12  to  24  hours  intervals.  The  displacement  vectors  were  generated  by 
tracking  perturbations  that  had  at  least  a  four  degree  difference  with  respect  to  the 
surrounding  sea  surface  temperature.  Vastano  and  Reid  proposed  that  these 
perturbations  (which  were  small  relative  to  the  baroclinic  radius  of  deformation)  were 
created  by  small-scale  wave  or  turbulent  processes  and  that  their  displacement 
represented  the  local  flow  regime.  Vastano  and  Reid  stated  that  Rossby-like  or 
internal  gravity  wave  motion  could  be  present,  but  that  the  wave  phase  speeds  would 
be  much  less  than  maximum  theoretical  values  and  therefore  assumed  that  these  would 
move  slower  than  typical  mesoscale  surface  currents.  An  additional  assumption  that 
Vastano  and  Reid  made  was  that  diffusion,  heat  exchange  with  the  atmosphere  and 
vertical  motions  were  much  less  than  the  exchange  of  heat  induced  by  advection  over 
the  time  period  analyzed.  The  generated  vectors  were  then  used  to  compute 
streamfunctions  which  allowed  the  nonuniform  vector  distributions  to  be  transformed 
into  a  uniform  vector  field.  Vastano  and  Reid  have  shown  that  vectors  generated  in 
this  manner  agree  quite  well  with  estimated  sea  surface  topography  in  the  area  oi~  the 
Oyashio  Front  northeast  of  Japan.  However,  the  primary  disadvantage  of  this  method 
is  that  it  requires  cloud-free  views  of  the  sea  surface. 

In  an  effort  to  move  away  from  the  subjectiveness  of  feature  tracking,  Emery  et 
al.  (19S6)  developed  a  method  that  was  objective  and  operator-independent  in 
computing  surface  advection  velocities  from  sequences  of  satellite  images.  In  this 
method  the  How  vectors  were  generated  by  identifying  the  end  points  of  the 
displacement  vectors  as  determined  by  the  maximum  cross  correlation  between  a  32  x 
32  pixel  search  window  in  the  second  image  and  a  22  x  22  pixel  template  in  the  first 
image.  Once  the  satellite  images  were  transformed  into  maps  of  SST  gradient,  the 
maximum  cross  correlation  (MCC)  technique  was  applied.  The  resulting  vectors  were 
spatially  coherent  and  exhibited  excellent  agreement  with  tracks  of  shallow-drogued 
drifters.  However,  near  coastal  boundaries  there  was  a  problem  in  that  some  of  the 
windows  contained  land  and  would  not  allow  consistent  representations  of  the 
advective  SST  changes  to  be  generated.  As  a  result,  all  such  velocity  vectors  were 
edited  out  of  the  final  vector  fields.  Another  problem  that  was  detected  with  the  MCC 
method  was  that  it  could  not  respond  to  rotational  or  deformational  changes  while  the 
feature-tracking  technique  could. 
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As  shown  above,  each  method  has  demonstrated  that  satellite-generated  velocity 
vectors  can  be  used  to  determine  surface  advection  velocities,  and  yet  each  method  also 
has  some  limitation.  To  what  extent  these  techniques  could  provide  accurate  ocean 
surface  velocities  in  the  coastal  region  is  unknown.  Therefore,  while  acknowledging  the 
subjectiveness  of  the  feature-tracking  technique,  this  method  was  chosen  because  no 
one  to  date  has  thoroughly  examined  what  ocean  motion  these  vectors  represent  in  the 
coastal  ocean  regime.  Additionally,  the  programming  required  to  develop  such  a 
method  was  straightforward  and  easily  implemented. 
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III.  TECHNIQUE 

For  this  research  AVIIRR  data  (in  the  form  of  ocean  flow  vectors)  was 
compared  against  in  situ  Doppler  current  measurements  and  hydrographic  data  which 
were  acquired  during  the  April  1981  phase  of  the  Coastal  Ocean  Dynamics  Experiment 
(CODE)  research  off  the  coast  of  California.  The  results  of  the  CODE  research  (Huyer 
and  Kosro.  1987)  were  used  as  ground  truth  for  all  subsequent  comparisons.  The 
satellite  data  was  acquired  at  Scripps  Institution  of  Oceanography  and  those  passes 
which  covered  the  area  of  interest  during  the  period  of  the  CODE  research  were  used 
to  create  the  ocean  flow  field.  The  April  1981  time  frame  was  chosen  after  viewing  the 
available  satellite  images  and  determining;  1)  which  images  had  the  least  amount  of 
cloud  contamination  in  the  CODE  research  region  and  2)  which  images  coincided  with 
one  of  the  periods  examined  by  Huyer  and  Kosro. 

A.       SATELLITE  IMAGERY 

1.  Description 

The  satellite  images  that  were  used  to  generate  the  ocean  flow  vectors  were 
512  X  512  pixel  portions  of  three  consecutive  XOAA-6  orbits.  The  data  were 
calibrated  in  aperture  brightness  temperature  and  resampled  onto  a  so-called  "equi- 
rectangular"  map  grid  (which  was  equivalent  to  creating  a  small  area  plotting  sheet 
using  a  "plane-sailing  approximation").  The  error  associated  with  this  assumption  will 
be  addressed  in  a  later  section.  The  images  were  loaded  onto  the  Naval  Postgraduate 
School  (NPS)  COMTAL  system  and  examined  to  determine  the  extent  and  location  of 
any  cloud  coverage  which  would  contaminate  the  infrared  signal  the  satellite  received. 
The  result  of  the  comparison  between  the  visible  (VIS)  and  the  infrared  (IR)  images  for 
the  first  pass  (04,27,81  16:15:00  GMT)  was  that  some  thin  high  clouds  were  detected 
near  the  western  edge  of  the  CODE  region  (Figure  3.1).  For  the  middle  (night)  pass 
(04  28  81  03:34:00  GMT)  some  low  stratus  along  the  coast  was  detected  by  noting  the 
changes  in  the  coastline  (Figure  3.2).  Additionally,  some  mid  and  high  clouds  were 
also  present,  and  these  limited  some  feature  tracking.  Comparison  of  the  VIS  and  IR 
images  for  the  third  pass  (04/28/81  15:52:00  GMT)  revealed  only  a  very  small  amount 
of  low  stratus  along  the  western  edge  of  the  research  area  (Figure  3.3). 
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Figure  3.1     XOAA-6  AVHRR  Infrared  image  04,  27  81  16:15:00  GMT. 
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Figure  3.2    NOAA-6  AVHRR  Infrared  image  04/28  81  03:34:00  GMT. 
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Fisure  3.3     NOAA-6  AVHRR  Infrared  image  04  28  81  15:52:00  GMT. 
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2.  Error  Analysis 

The  plane-sailing  approximation  can  he  described  as  centering  an  equi- 
rectangular  plane  tangent  to  the  earth  at  a  known  point  in  the  image(s).  For  all  three 
passes  the  tangent  point  was  chosen  to  be  Point  Reyes,  with  a  known  position  of 
3S°N,  123°\V.  Once  the  center  point  of  the  plane  had  been  selected,  the  image  was 
then  projected  onto  the  equi-rectangular  grid  of  the  plane.  To  determine  the  error 
associated  with  this  approximation,  twelve  landmarks  were  selected  and  their  line  and 
pixel  locations  determined.  The  line  and  pixel  locations  were  then  converted  to  image 
latitudes  and  longitudes  based  on  location  of  Point  Reyes  in  each  image.  These  image 
latitudes  and  longitudes  were  then  compared  against  the  true  latitudes  and  longitudes 
which  were  determined  using  U.S.  Department  of  Commerce  Sectional  Aeronautical 
Charts  for  the  region  of  interest.  The  differences  between  the  image  and  the  true 
latitudes  and  longitudes  revealed  two  important  pieces  of  information. 

First,  all  the  images  needed  to  be  rotated  approximately  8  degrees  to  the  left 
(counter-clockwise)  in  order  to  align  each  image  with  true  North.  Unfortunately,  this 
correction  could  not  be  made  on  the  COMTAL  due  to  software  limitations  (the 
COMTAL  would  only  allow  image  rotations  in  90  degree  increments).  However,  this 
correction  was  applied  through  hand  calculation  and  all  the  coastal  landmarks  then 
plotted  within  2  pixels  of  their  true  location.  This  residual  error  (approximately  7  -  10 
cm  sec)  can  be  traced  to  the  geometry  between  the  satellite,  the  view  angle,  and  the 
earth  location  of  the  pixel(s).  While  rotating  the  images  eliminated  some  of  the 
navigational  error  due  to  applying  the  plane-sailing  approximation,  the  remaining  error 
of  7  -  10  cm  sec  cannot  be  corrected  without  solving  for  the  associated  rotational  effect 
on  the  different  spherical  triangles  involved  in  the  problem.  Unfortunately,  these 
calculations  could  not  be  made  because  a  program,  (now  being  developed  at  the  Naval 
Postgraduate  School)  which  navigates  an  image  by  obtaining  a  best  fit  for  several 
landmarks  via  changing  the  satellite's  roll,  pitch,  and  yaw,  was  not  available  at  the  time 
this  study's  AVHRR  data  were  processed.  (For  more  information  on  this  program  see 
the  recommendations  section  of  Chapter  V.)  Therefore,  it  was  necessary  to  reduce  the 
image  navigation  error  via  the  feature  tracking  program  (Oceantrak)  as  explained  in 
the  next  section. 

The  second  important  piece  of  information  that  was  expected  but  could  not  be 
determined  precisely  was  the  East;  West  error  due  to  using  an  equi-rectangular  grid 
when  lines   of  longitude   on  the  earth   converge  as   the  distance   from  the  equator 
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increases.  Since  Point  Reyes  was  selected  as  the  center  point  in  processing  these 
images,  the  assumed  value  of  1.1  km  per  pixel  would  be  valid  for  any  pixel  along  the 
38th  parallel  of  latitude.  North  of  this  latitude,  the  "width''  of  each  pixel  should  get 
smaller  by  approximately  15  meters  per  degree  of  latitude.  Similarly,  south  of  3S°\ 
each  pixel  should  get  larger  by  approximately  the  same  amount  (15  meters)  per  degree 
of  latitude.  However,  in  the  region  of  this  study  (Point  Reyes  to  Point  Arena)  the 
difference  in  pixel  size  would  only  change  the  computed  velocities  by  a  maximum  of 
0.037  cm/sec. 

B.       OCEANTRAK 

1.  Description 

A  modified  cloud  tracking  program  (Oceantrak)  was  used  to  generate  the  flow 
vectors.  Oceantrak  required  the  user  to  input  a  "crossing"  time  for  each  pass  so  that  a 
scaling  factor  cm/ps  (centimeters  per  pixel  per  second)  could  be  computed.  The 
resulting  factor  was  then  used  to  calculate  a  velocity  based  on  the  displacement  of  a 
submesoscale  feature  seen  in  both  images. 

The  displacement  vectors  were  created  by  "marking"  a  feature  with  a  movable 
cursor  on  one  image,  and  then  "toggling"  to  the  second  image  and  "marking"  the  same 
feature  again.  Because  a  water  vapor  correction  was  not  made  to  any  of  the  three 
images,  the  gray  shade  of  most  "marked"  features  changed  between  images.  To  avoid 
"marking"  different  features,  unique  bends,  bumps,  and  other  small  scale  perturbations 
were  chosen  to  be  tracked.  Offshore  these  types  of  features  were  quite  easily  tracked, 
unless  cloud  cover  moved  into  the  area.  But,  as  the  coast  was  approached,  it  became 
increasingly  difficult  to  locate  the  same  feature  from  one  pass  to  another,  and  therefore 
the  number  of  vectors  generated  decreased  as  the  coast  was  approached.  Additionally, 
features  that  did  not  appear  to  move  were  tracked  by  "marking"  the  center  of  the 
feature  in  both  images.  However,  only  a  few  of  these  types  of  features  were  present  in 
the  images  that  were  examined. 

Oceantrak  was  also  used  to  reduce  the  extraneous  movement  of  features  and 
landmarks  due  to  an  initial  lack  of  image  to  image  coregistration.  Because  successive 
satellite  paths  are  not  the  same,  the  images  must  be  aligned  (coregistered)  so  that 
landmarks  will  coincide  spatially  on  the  screen.  Since  neither  of  the  sets  of  passes  were 
coregistered.  the  middle  pass  was  chosen  to  serve  as  ground  truth  for  registering  the 
other  two  passes.    The  images  were  registered  by  first  determining  the  movement  of 
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several  land  features  and  then  attempting  to  minimize  the  movement  in  the  area  of  the 
CODE  data.  By  using  Oceantrak  to  coregister  the  images  most  of  the  error  from  the 
plane-sailing  approximation  was  eliminated.  Because  the  images  could  only  be  moved 
by  integer  pixel  amounts,  there  was  some  slight  residual  movement  (error)  between 
images. 

2.  Error  Analysis 

For  this  experiment  the  satellite's  "crossing''  time  over  San  Francisco  Bay  (to 
the  nearest  minute)  was  determined.  The  maximum  error  resulting  from  determining 
the  "crossing"  time  to  the  nearest  minute  was  0.004  cm  p  s  (less  than  0.2  percent). 
While  the  error  associated  with  the  "crossing"  time  was  quite  small,  the  error  due  to 
imperfect  coregistration  of  the  images  ranged  from  2.7  cm  sec  from  the  East  for  the 
first  set  of  vectors  to  3.5  cm  sec  from  the  Northeast  for  the  second  set.  These  last  two 
errors  are  due  to  the  limitation  created  by  only  being  able  to  move  the  images  by  an 
integer  pixel  amount. 

C.       CODE  DATA 

1.  Description 

The  CODE  Doppler  profile  data  was  acquired  through  the  continuous  use  of 
an  Ametek-Straza  Doppler  acoustic  log  (DAL).  The  horizontal  velocity  vector  profiles 
spanned  the  15-150  m  depth  range  and  the  data  points  were  vertically  separated  by  6.5 
m.  The  currents  were  estimated  by  filtering  the  combined  DAL  and  navigation 
(LORAN-C)  data  over  30  minute  intervals.  The  CODE  survey  grid  was  designed  to  be 
small  enough  to  be  completed  in  2-3  days,  and  yet  large  enough  to  cover 
approximately  100  km  of  coastline.  Since  the  time  required  to  complete  a  survey  was 
substantially  less  than  the  typical  duration  (7  to  10  days)  of  a  strong  wind-driven 
upwelling  event,  each  of  the  mesoscale  surveys  was  considered  "synoptic",  or  at  least 
quasi-synoptic. 

2.  Error  Analysis 

For  this  study  the  Doppler  current  vectors  were  obtained  from  the  Huyer  and 
Kosro  (1987)  paper.  The  errors  associated  with  these  vectors  were  determined  to  be 
2-3  cm;  sec  in  the  alongshore  component  and  4-5  cm; sec  in  the  cross-shore  component 
(Kosro,  1985).  That  portion  of  the  signal  due  to  inertial  motion  or  internal  tides  is  not 
quantified  in  Huyer  and  Kosro  (1987).  The  potential  influence  this  submesoscale 
variability  might  introduce  was  addressed  in  another  survey  (Kosro  and  Huyer,  1986). 
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After  conducting  surveys  on  two  jets  oil  the  coast  of  northern  California, 
Kosro  and  Huyer  (19S6)  concluded  that  the  offshore  jets  were  likely  continuations  of 
the  associated  nearshore  jets.  To  determine  the  submesoscale  variability  (caused  by 
internal  tides  and  inertial  motion)  in  the  DAL  current  and  hydrographic  measurements, 
Kosro  and  Huyer  separated  the  (expected)  geostrophic  mesoscale  flow  from  the  smaller 
scale  and  more  rapidly  fluctuating  components.  This  was  accomplished  by  using 
objective  analysis  to  estimate  the  geostrophic  current  at  each  DAL  vector  location 
beyond  the  500  meter  isobath,  and  then  both  sets  of  vectors  were  smoothed  using  a 
Gaussian  filter.  Once  the  small  scale  fluctuations  had  been  removed,  Kosro  and  Huyer 
found  that  the  correspondence  between  the  geostrophic  and  the  DAL  currents  was 
"striking".  Assuming  that  the  above  jets  were  truly  extentions  of  the  associated 
nearshore  jets,  then  ignoring  internal  and  inertial  motion  in  the  nearshore  region 
should  not  introduce  any  large  errors.  For  this  reason  I  have  not  included  them  in  the 
error  budget  table. 

D.       TOTAL  ERROR  BUDGET 

The  largest  source  of  error  in  this  study  was  the  use  of  an  inaccurate  navigation 
routine  in  navigating  the  images.  This  method  of  navigation  and  the  resulting 
imperfect  coregistration  of  the  two  images  introduced  errors  of  up  to  2.7  cm  sec  in  the 
computation  of  all  velocity  vectors  generated  using  the  first  and  second  images. 
However,  the  resulting  imperfect  coregistration  from  using  this  method  of  navigation 
introduced  errors  ranging  as  large  as  3.5  cm/sec  in  the  generation  of  the  second  set  of 
velocity  vectors.  While  these  errors  are  large  (in  comparison  to  the  error  in  the  pixel 
size  due  to  the  converging/diverging  line  of  longitude),  a  qualitative  comparison 
between  the  Doppler  data  and  the  satellite-derived  vectors  is  still  be  possible  since  the 
error  in  the  DAL  data  is  as  large  or  larger.  See  Table  1  for  a  breakdown  of  the 
different  error  components. 
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TABLE  1 

TOTAL  ERROR  BUDGET 

ERROR 
SOURCE 

NORTH  -  SOUTH 
ERROR  ( cm/ sec ) 

EAST  -  WEST 
ERROR  ( cm/sec) 

Plane  Sailing 
Approximation 

. 00023-. 00025 

0.  034  -  0.  037 

Crossing  Time 

0.  07  -  0. 08 

0.  07  -  0.  08 

Imoerf ect 

Cor egi strati on 

Pass  1-2 

0.  00 

2.  70 

Imperfect 

Co registration 

Pass  2-3 

2.  48 

2.  48 

Dopp  ier 
Acoustic  Log 
(  ?Cosro,  198 o  ) 

2.0  -  3.0 

4.  0  -  5.  0 
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IV.  RESULTS 

As  expected  from  California  Current  climatology,  the  general  flow  in  the  CODE 
region  was  towards  the  south  during  the  April  19S1  survey.  The  dynamic  topography 
as  determined  from  the  data  acquired  during  the  April  survey  also  agreed  (in  general 
terms)  with  the  climatological  mean  dynamic  topography  for  the  month  of  April 
(Hickey,  1979).  The  winds  during  the  April  survey  were  strongly  favorable  for 
upwelling  (Huyer  and  Kosro.  19S7). 

A.       CODE  SURVEY 

All   of  the   following  descriptions  of  hydrographic  and   DAL  features  are  my 
interpretations  of  data  presented  in  Huyer  and  Kosro  (19S7). 
1.  Hydrographic  Features 
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Figure  4.1     Maps  of  in  situ  sea  surface  temperature  and  temperature 
at  2U  m  (from  Huyer  and  Kosro.  1987). 

Plots  of  the  April  survey  SST  showed  the  coldest  surface  water  adjacent  to  the 
coast,  and  a  strong  offshore  temperature  gradient,  both  of  which  are  consistent  with 
the  strong  upwelling  favorable  winds  (Figure  4.1).  Maps  of  the  surface  salinity  also 
reinforce  the  fact  that  upwelling  was  occurring  during  the  study  period,  in  that  the 
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more  saline  water  was  located  adjacent  to  the  coast  (Figure  4.2).  Both  the  temperature 
and  salinity  distributions  at  20  meters  also  verify  that  upwelling  was  taking  place 
during  the  April  survey.  The  map  of  geopotential  anomaly  of  the  sea  surface  relative 
to  100  dbar  was  very  similar  to  the  map  of  the  near-surface  salinity,  with  the  highest 
values  being  associated  with  the  low-salinity  water  offshore  and  the  lowest  values 
coinciding  with  the  high-salinity  water  near  the  coast  (Figure  4.3).  Iluyer  and  Kosro 
propose  that  southward  advection  of  subarctic  water  would  explain  the  persistently  low 
surface  salinities  offshore,  while  northward  advection  of  local  runoff  through  a  narrow 
and  sporadic  countercurrent  along  the  coast  would  explain  the  isolated  pockets  of  low 
salinity  water  over  the  shelf.  Huyer  and  Kosro  also  state  that  the  circulation  pattern 
they  observed  during  the  survey  was  much  more  complex  than  they  had  expected.  In 
fact,  on  the  northern  edge  of  the  survey  area,  all  of  the  above  maps  indicate  the 
presence  of  a  warm,  low-salinity,  anticyclonically  rotating  eddy.  A  vertical  phase  shift 
of  this  eddy  can  be  seen  in  the  temperature  maps  (see  Figure  4.1).  The  surface  position 
of  the  eddy  as  shown  by  the  SST  map  agrees  well  with  the  satellite  imagery.  The  20 
meter  position  of  the  eddy  (which  is  southwest  of  the  surface  position)  as  shown  by  the 
temperature  plot  at  that  level  agrees  with  the  plot  of  the  DAL  data  at  the  20  meter 
level.  In  the  southwest  region  of  the  study  area  there  is  little  indication  (due  to  the 
spacing  of  the  hydrographic  stations)  of  a  second  warm,  low-salinity,  anticyclonically 
rotating  eddy  present  in  the  above  maps.  Unfortunately,  a  comparison  of  the  surface 
and  20  meter  positions  of  this  eddy  cannot  be  made  due  to  the  coarse  resolution  of  the 
hydrographic  data  in  the  area  of  interest. 
2.  DAL  Features 

The  20  meter  map  of  the  DAL  current  vectors  (Figure  4.4)  for  the  April 
survey  shows  the  southern  eddy  quite  well,  but  does  not  give  as  good  a  representation 
of  the  eddy  to  the  north.  The  DAL  map  at  20  meters  also  shows  a  general  southward 
flow  in  the  region. 

At  80  and  150  meters,  the  DAL  current  velocity  maps  (Figure  4.4)  are  very 
similar  to  the  20  meter  map.  However,  as  expected,  the  number  of  near  shore  velocity 
vectors  decreases  as  the  water  shoals.  Offshore,  both  of  the  eddies  are  represented  to 
the  same  degree  as  they  are  in  the  20  meter  map  and  the  general  southward  How  of  the 
region  is  also  still  evident.  The  only  real  differences  between  the  20,  80,  and  150  meter 
maps  are  (as  expected)  that  the  magnitude  of  all  the  velocities  decreases  as  the  depth 
increases. 
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Figure  4.2     Maps  of  in  situ  surface  salinity  and  salinitv  at  20  m 
(from  Huyer  and  Kosro,  F9S7). 


Fisure  4.3     Map  of  the  seopotential  anomalv  (m~  s") 
of  the  sea  "surface  relative  to  1  Rf dbar  {from  Huyer  and  Kosro,  19S7). 
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Fieure  4.4     Mars  of  DAL  current  vectors  at  depths  of  20.  SO.  and 
15u  m  (irom  Huyer  and  Kosro.  1987). 

B.       SATELLITE  DATA 

In  the  study  region,  the  one  feature  that  stands  out  quite  clearly  is  the  large 
warm  eddy  along  the  northern  edge  (Figure  4.5).  Although  this  eddy  was  partially 
covered  by  clouds  in  the  middle  (or  night)  pass,  the  majority  of  the  feature  was 
detectable,  particularly  the  boundaries.  Unfortunately,  the  second  eddy  to  the  south 
was  not  as  detectable.  In  fact  the  southern  edge  of  this  feature  was  quite  diffuse  and 
as  a  result,  a  boundary  could  not  be  determined  from  the  SST  imagery.  This 
homogeneity  of  the  coastal  water  mass  was  also  a  problem  in  generating  near-shore 
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velocity  vectors.  The  near-shore  problem  was  compounded  by  coastal  fog  that 
developed  during  the  night  and  severely  limited  use  of  the  middle  pass  in  generating 
velocity  vectors  near  the  coast. 

Outside  the  study  area  (off  Monterey  Bay)  there  were  several  distinct  SST 
features  that  were  easily  tracked  and  which  appeared  in  all  three  images  used  in  this 
study.  Velocity  vectors  were  generated  in  this  area  in  the  hope  that  another  in  situ 
data  set  could  be  found  to  serve  as  a  ground  truth.  Unfortunately,  I  could  not  find 
such  a  data  set.  These  other  vectors  were  left  on  the  images  so  that  the  reader  could 
see  how  fast  some  of  the  features  moved,  and  more  importantly  how  active  the  coastal 
region  is  on  a  mesoscale  basis  (Figures  4.5  and  4.6). 

The  SST  features  in  the  study  region,  as  shown  by  the  second  set  of  satellite 
images,  became  quite  diffuse.  This  is  due  in  part  to  the  fog  which  developed  during  the 
night  and  obscured  many  of  the  nearshore  features.  Another  factor  is  that  a  large 
cloud  mass  can  be  seen  moving  into  the  study  region  from  the  northwest.  This  cloud 
mass  and  its  associated  water  vapor  also  contaminated  a  portion  of  the  satellite  sensed 
data  for  the  region.  Because  of  these  factors,  approximately  half  as  many  velocity 
vectors  were  generated  using  the  second  set  of  images  as  were  generated  using  the  first 
set  (Figure  4.6). 

C.       COMPARISON  OF  CODE  AND  SATELLITE  VECTORS 

Because  only  ten  vectors  were  generated  from  the  second  set  of  images,  the 
following  comparison  will  concentrate  on  those  vectors  generated  from  the  first  set  of 
satellite  images,  where  nineteen  vectors  were  generated.  Figure  4.7  shows  the  first  set 
of  satellite  vectors  (large  arrowheads)  superimposed  on  the  20  meter  DAL 
measurements. 

1.  Agreement 

Overall,  the  first  set  of  satellite  velocity  vectors  that  were  generated  in  the 
study  area  agree  quite  well  (in  direction  of  How)  with  the  DAL  vectors  (Figure  4.7). 
Where  comparisons  could  be  made  (both  types  of  data  located  in  approximately  the 
same  position  and  the  directions  of  flow  agreed),  the  DAL  measurements  were 
approximately  1.5  times  larger  than  those  of  the  associated  satellite  derived  velocity 
vectors.  Eight  out  of  eleven  vectors  agreed  with  DAL  measurements  where  the  two 
data  sources  coincided.  An  additional  eight  satellite  vectors  were  generated  at 
locations  in  the  study  area,  but  outside  the  domain  of  the  DAL  measurements  so  that 
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Figure  4.5     Satellite  vectors  generated  from  passes  1  and  2. 
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Figure  4.6     Satellite  vectors  generated  from  passes  2  and  3. 
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Figure  4.7     Map  of  20  m  DAL  and  satellite-derived  velocity  vectors. 
Larec  arrowheads  are  satellite  venerated  vectors, " 
small  arrowheads  are  DAL  measurements. 

comparisons  could  not  be  made.  The  anticyclonic  rotation  of  the  northern  eddy  is 
verified  by  the  satellite  vectors,  as  well  as  the  general  southward  How  of  the  region. 
The  satellite  vectors  also  show  some  support  for  the  southern  eddy  and  show  rotation 
in  the  correct  direction  (anticyclonically). 

Although  Huyer  and  Kosro  (19S7)  proposed  a  sporadic  countercurrent  along 
the  coast  could  explain  the  isolated  pockets  of  low  salinity  water  over  the  shelf,  the 
DAL  measurements  do  not  reveal  any  countercurrent.    Similarlv,  the  satellite-derived 
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vectors  do  not  indicate  the  presence  of  any  countercurrent  along  the  coast.  In  fact. 
nearshore  both  the  DAL  measurements  and  the  satellite  generated  velocity  vectors 
agree  quite  well. 

2.  Disagreement 

There  are  three  satellite-derived  velocity  vectors  which  disagree  with  the  DAL 
measurements  (Figure  4.7).  The  first  is  due  to  frontal  waves  propagating  southward 
along  the  sharp  boundary  between  the  cold  upwelled  water  and  the  warmer  offshore 
water.  These  waves  exist  along  the  northern  edge  of  the  study  region  (where  they  are 
180  degrees  out  of  phase  with  the  northward  DAL  data).  One  of  the  resulting  velocity 
vectors  (3  cm  sec,  which  is  near  the  noise  level  of  the  technique)  is  present  in  Figure 
4.7  This  vector  is  located  northwest  oi"  Pt.  Arena,  over  the  1000  meter  isobath.  The 
frontal  waves  also  follow  the  western  boundary  of  the  northern  eddy  (where  two 
additional  waves  were  tracked).  When  these  features  were  originally  tracked,  they  were 
so  clear  in  the  imagery  that  the  eddy  appeared  to  be  rotating  cyclonically.  It  was  only 
after  determining  that  the  features  were  not  part  of  the  eddy  but  rather  traveling  along 
the  western  boundary  that  separated  the  warm  offshore  water  from  the  cold  upwelled 
water  that  I  realized  they  had  to  be  frontal  waves. 

The  second  disagreement  between  the  satellite  vectors  and  the  DAL 
measurements  occurred  along  the  northern  edge  of  the  southern  eddy.  In  attempting 
to  track  the  diffuse  edge  of  the  eddy,  one  northward  vector  was  generated  which  was 
approximately  180  degrees  out  of  phase  with  the  southward  DAL  vector.  This  vector 
is  located  southwest  of  Pt.  Arena,  again  over  the  1000  meter  isobath  and  is  labelled  20 
cms.  This  northward  shift  in  the  location  of  the  eddy  is  most  likely  due  to 
atmospheric  contamination  from  a  small  group  of  clouds  that  moved  into  area  during 
the  night. 

There  is  a  third  vector  that  disagrees  with  the  DAL  measurements,  but  the 
velocities  of  both  vectors  are  small  (4  -  5  cms).  This  satellite  vector  is  located  over  the 
150  meter  isobath,  practically  due  west  of  Pt.  Arena  and  is  labelled  4  cm/s.  The  vector 
was  generated  by  tracking  what  appeared  to  be  a  perturbation  on  the  edge  of  the  warm 
eddy.  While  it  is  possible  that  the  perturbation  was  an  instability  on  the  eddy,  clouds 
had  moved  into  this  particular  area  between  the  two  passes  used  to  generate  the  vector. 
It  is  possible  that  the  perturbation  is  atmospheric  in  nature,  not  oceanic. 


34 


3.  Other  Physical  Processes 

Because  the  winds  during  the  April  1981  survey  were  upwelling  favorable,  one 
might  expect  significant  ofTshore  directed  velocity  vectors  associated  with  Ekman 
transport.  An  estimate  of  this  contribution  to  the  satellite  derived  motion  can  be 
obtained  from  simple  analytic  solutions  (see  Csanady,  1982).  Csanady's  solutions 
provide  cross-shelf  velocities  which  are  zero  at  the  coast  and  increase  exponentially 
ofTshore.  The  maximum  cross-shore  response  occurs  at  a  distance  offshore  equal  to 
the  external  Rossby  radius  of  deformation  (350  -  500  kilometers).  However,  Kelly 
(1983)  has  shown  that  coastal  topography  plays  a  significant  role  in  the  wind  forced 
response  of  the  circulation  in  the  CODE  region.  The  wind  driven  ocean  response  is 
thus  spatially  variable.  The  satellite  vector  fields  (Figures  4.5  and  4.6)  show  no 
obvious  uniform  offshore  velocity  trend.  Based  on  Csanady's  analytic  solutions,  strong 
ofTshore  velocities  are  suspected  to  occur  further  ofTshore. 

The  good  agreement  between  the  satellite-derived  vectors  and  the  subsurface 
geostrophic  flow  indicates  that  significant  contributions  by  inertial  motion  are  not 
present  in  the  satellite  vectors.  This  is  also  in  agreement  with  the  results  of  Kosro  and 
Huyer  (1986),  as  discussed  above. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

Satellite-derived  ocean  surface  velocities  have  been  compared  with  in  situ 
observations.  In  general,  the  satellite-derived  velocities  agree  with  the  hydrographic 
surveys  and  the  Doppler  acoustic  log  data.  This  implies  that  the  vectors  during  this 
time  period  are  representative  of  the  subsurface  geostrophic  How.  Additionally,  the 
satellite  vectors  do  not  show  any  signs  of  any  inertial  motion  or  Ekman  forcing. 

The  satellite-derived  velocity  vectors  agreed  with  the  hydrographic  data  that  was 
acquired  during  the  April  1981  CODE  survey.  But,  while  the  in-situ  hydrographic  data 
could  not  resolve  the  mesoscale  jet  and  eddies  located  in  the  study  region,  the  satellite 
generated  vectors  did  resolve  these  features.  The  satellite  vectors  were  found  to  be 
representative  of  the  subsurface  geostrophic  flow,  due  to  the  baroclinic  nature  of  the 
study  region.  While  this  technique  can  generate  velocity  vectors  which  are 
representative  of  the  geostrophic  How,  it  is  important  to  note  that  the  success  of  this 
technique  depends  on  the  water  vapor  content  of  the  atmosphere.  In  particular,  the 
effects  of  fog  (prevalent  during  during  upwelling  conditions)  may  preclude  the  use  of 
this  technique. 

Like  the  offshore  study  by  Flament  et  al.  (1985),  DAL  measurements  were 
approximately  1.5  times  larger  than  the  satellite-derived  velocity  vectors  (when  they 
were  co-located).  However,  the  number  and  concentration  of  the  satellite  vectors  were 
substantially  smaller  than  the  DAL  data.  This  sparseness  of  data  made  point  by  point 
comparisons  of  the  DAL  versus  satellite  data  difficult. 

B.  POTENTIAL  USES 

There  are  many  potential  uses  for  this  technique,  both  as  currently  demonstrated 
and  with  some  minor  improvements.  For  example,  if  the  concentration  of  velocity 
vectors  could  be  increased,  via  either  feature  tracking  from  color  images  or  by  applying 
an  atmospheric  correction  (so  that  isotherm  displacements  might  be  tracked),  the 
technique  could  definitely  provide  hypothetical  input  data  for  global  ocean  circulation 
modeling  and  possibly  enough  data  to  initialize  a  regional  ocean  circulation  model.  In 
its  current  operating  condition,  the  technique  could  be  of  assistance  in  SAR  (Search 
And  Rescue)  operations  by  allowing  the  user  to  visualize  ocean  surface  flows  (direction 
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and  speed)  in  any  area  of  interest.  This  information  could  then  be  used  to  organize  the 
search  to  take  maximum  advantage  of  assets  (aircraft  and  or  ships)  available  in 
locating  the  victim  in  distress.  Information  on  current  variability  every  6-12  hours 
would  also  be  quite  valuable  in  tracking  pollution  (oil  spills)  and  would  assist  in 
predicting  if,  when  and  possibly  where,  such  pollution  might  come  ashore.  The  lead 
time  this  information  would  provide  could  then  be  used  to  notify  local  residents  and 
position  clean-up  equipment  in  those  areas  where  it  would  do  the  most  good.  This 
technique  could  also  be  used  in  ASW  ( Anti-Submarine  Warfare)  by  allowing  features 
of  interest  (fronts  and  eddies)  to  be  tracked  and  their  location  and  movement  briefed  to 
those  forces  in  the  associated  region.  On  a  more  scientific  note,  this  same  information 
would  be  very  useful  in  assisting  shipboard  hydrographic  surveys  o[  both  fronts  and 
eddies  by  providing  the  ship(s)  with  close  to  real-time  (shortly  after  the  satellite 
image(s)  have  been  processed)  location  and  movement  of  such  features  and  by  allowing 
gaps  in  the  in  situ  data  to  be  filled  via  the  information  obtained  from  the  velocity 
vectors.  The  above  are  just  a  few  of  the  potential  uses  of  this  and  similar  techniques; 
undoubtly  more  uses  and  more  accurate  information  will  be  found  as  research  in  this 
field  continues. 

C.       RECOMMENDATIONS 

In  order  to  obtain  more  accurate  information  from  this  technique,  the  satellite 
images  that  are  used  must  be  navigated  to  the  highest  degree  possible.  A  new  AVHRR 
navigation  program  is  presently  being  developed  at  the  Naval  Postgraduate  School  (Dr. 
James  L.  Mueller  personal  communication,  1987).  This  program  minimizes  the  error  in 
the  location  of  several  user  selected  landmarks  by  adjusting  the  satellite's  roll,  pitch, 
and  yaw.  Once  the  overall  error  has  been  minimized,  the  image  is  resampled  and 
remapped  with  a  latitude  and  longitude  computed  for  every  sixteenth  line  and  pixel. 
This  approach  to  navigating  satellite  images,  together  with  image  resampling  onto  a 
conformal  map  projection  (i.e.  transverse  Mercator)  would  significantly  reduce  the 
distortion  which  was  the  dominant  source  of  error  in  determining  displacement  vectors 
in  this  study.  Considering  the  above,  it  is  recommended  that  if  additional  research  is 
conducted  on  this  technique  that  the  AVHRR  navigation  program  be  used  when  it 
becomes  available. 

An  additional  recommendation  for  further  research  with  this  technique  is  that 
color  enhanced  images  might  provide  a  denser  concentration  of  velocity  vectors  by 
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possibly  allowing  smaller  temperature  differences  to  be  tracked.  However,  to  color 
enhance  the  satellite  images  properly,  an  atmospheric  correction  will  have  to  be  made 
to  allow  an  accurate  determination  of  the  SST.  This  correction  would  be  necessary  to 
avoid  color  shifts  between  images  that  are  only  the  result  of  an  increase  or  decrease  in 
water  vapor  in  the  atmosphere.  Since  this  correction  requires  two  IR  channels,  only 
satellite  data  from  NOAA-7  and  odd  numbered  follow-on  satellites  can  be  used. 
Fortunately,  there  are  additional  surveys  (during  other  periods)  contained  in  the  Huyer 
and  Kosro  (1987)  report  to  allow  a  comparison  to  be  conducted  using  color  tracking 
(NOAA-7  images)  and  in-situ  data. 

There  are  two  potentially  fruitful  areas  of  study  which  only  require  one  IR 
channel's  worth  of  data.  The  first  recommendation  is  to  investigate  a  no  wind  or  low 
wind  situation  to  determine  the  extent  of  influence  the  wind  has  or  does  not  have  on 
the  concentration  and  magnitude  of  the  satellite-derived  velocity  vectors.  Again,  the 
Huyer  and  Kosro  (1987)  report  contains  information  acquired  during  three  surveys  in 
which  the  winds  were  quite  weak.  This  would  allow  additional  satellite  generated 
velocity  vectors  to  be  compared  with  in  situ  hydrographic  data.  The  second 
recommendation  is  to  investigate  an  offshore  region  in  an  effort  to  determine  if 
satellite-derived  velocity  vectors  agree  or  disagree  with  in  situ  hydrographic  data. 
Although  a  study  similar  to  the  above  recommendation  has  been  mentioned  (Flament 
et  al.,  19S5),  a  detailed  comparison  between  in  situ  data  and  satellite  generated  velocity 
vectors  in  a  coastal  transition  region  has  not  been  published  as  far  as  this  author 
knows. 

Finally,  a  word  of  advice,  to  maintain  an  oceanographic  frame  of  mind,  avoid 
generating  vectors  just  for  the  sake  of  generating  vectors  and  determine  if  each  vector 
is  consistent  with  the  local  thermal  contrast  in  SST  and  geostrophy.  By  maintaining 
the  above  thought  process,  the  user  will  avoid  generating  false  vectors  associated  with 
large  scale  wave  propagation.  It  should  be  noted  that  the  above  reasoning  will  provide 
good  guidance  in  generating  velocity  vectors  in  most  cases.  When  salinity  becomes  as 
important  as  temperature  in  determining  the  density,  the  above  reasoning  would  no 
longer  apply.  This  type  of  situation,  in  which  salinity  appeared  to  be  at  least  as 
important  as  temperature  in  determining  surface  density,  was  encountered  by  Kosro 
and  Huyer  (1986)  in  an  offshore  study.  Therefore,  any  potential  users  of  this  or  other 
vector  generating  system  should  have  a  good  working  knowledge  about  the 
oceanographic  conditions  that  should  be  expected  in  the  region  where  they  are 
applying  this  technique. 
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